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A non-interacting electron gas on a one-dimensional ring is considered at finite temperatures. 
The localized spin is embedded at some point on the ring and it is assumed that the interaction 
between this spin and the electrons is the exchange interaction being the basis of the Ruderman- 
Kittel-Kasuya-Yosida indirect exchange effect. When the number of electrons is large enough, it 
turns out that any small but finite interaction radius value can always produce an essential change 
of the spin density oscillations in comparison with the zero interaction radius traditionally used to 
model the Ruderman-Kittel-Kasuya-Yosida effect. 



PACS numbers: 75.20.Hr, 75.75.+a, 71.10.-w 



I. INTRODUCTION 

The Ruderman-Kittel-Kasuya-Yosida (RKKY) effect 
is a phenomenon which plays an important role in a 
formation of magnetic structures in different systems. 
The essence of the effect consists in an indirect inter- 
action between localized spins placed in a Fermi sea of 
non-interacting conduction electrons. The interaction is 
called indirect because the spins feel the presence of each 
other through the electrons surrounding them: a local- 
ized spin interacting with the electrons induces in the 
electron gas spin density oscillations which then make 
impact on another localized spin. The interaction under- 
lying the RKKY effect, that is the interaction between 
the localized spin and the electron gas, can be of different 
nature. It can be either the hyperfme interaction between 
a localized nuclear spin and conduction electrons^ or the 
exchange interaction between the conduction and local- 
ized electrons^. The latter case is realized for example 
in alloys with transition metal ions where the conduction 
s-electrons and the localized c£-electrons of an ion interact 
through the s — d exchange interaction. 

To model the RKKY effect it is traditionally assumed 
that the interaction between a localized spin and conduc- 
tion electrons is local in the real space. This is modeled 
by Dirac's delta function^. In the case of the hyperfme 
interaction this model looks quite plausible because the 
size of a nucleus, being of order 10 -6 nm (see Ref. @), is 
small enough. The RKKY indirect exchange effect based 
on the hyperfme interaction was studied in the scientific 
literature. In Ref. the RKKY interaction between nu- 
clear spins embedded in a mesoscopic ring and a finite 
length quantum wire was investigated in the presence 
of a magnetic field. The indirect nuclear spin interac- 
tion was found to depend on the nuclear spin positions, 
number of the conduction electrons, magnetic field and 
system's geometry. The influence of electron-electron in- 
teractions and electron spin correlations on the RKKY 
interaction between two nuclear spins was considered by 
Semiromi et al& The nuclear spins were embedded in a 
mesoscopic metallic ring. It was numerically shown that 
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FIG. 1: (Color online) A non-interacting electron gas on a ID 
ring of radius R. The localized spin is centered around the 
ring point with the polar angle cj> = 0. The localization radius 
of the spin is ro which is also estimated to be the radius of 
the exchange interaction in the RKKY effect. 



the electron-electron interactions and electron spin cor- 
relations can essentially change the RKKY interaction 
dependence on the magnetic flux. 

However, in the case of the exchange interaction the 
ionic spins are much less localized. For example, the ionic 
radius, which we will denote by ro, for the /-shell metal 



ions Er 3+ and Nd 3H 



is r 



0.096 nm and r = 0.108 



nm (see Ref. l9j) , respectively. The value of ro gives an ef- 
fective radius of the exchange interaction. Moreover, one 
can easily conceive a situation where artificial objects like 
quantum dots with non-zero total spin are embedded in 
a Fermi sea of conduction electrons. The interaction be- 
tween the total spin of such objects and the electrons 
is similar to the exchange interaction and produces the 
RKKY interaction between the total spins of those arti- 
ficial objects. The value of ro can thus be technologically 
varied. Systems with quantum dots interacting through 
the RKKY effect were already investigated in a number 
of scientific papers. The RKKY effect between two quan- 
tum dots embedded in an Aharonov-Bohm ring was in- 
vestigated by Utsumi et al^ as a function of a magnetic 
flux through the ring. The quantum dots contained odd 
numbers of electrons. The interaction between the total 
quantum dots' spins and the conduction electrons in the 
ring was described by a tunneling Hamiltonian. In this 
tunneling Hamiltonian the coupling constant was mod- 
eled by the delta function that is the interaction radius 
was ro = 0. In Ref. 11 two localized spins in one-, two- 
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FIG. 2: (Color online) The function n s (0,ro) describing the 
spin density oscillations on a mesoscopic ring for the con- 
ventional model with ro = 0. The number of electrons is 
iVe] = 24. 



and three-dimensional electron gases were considered. 
Decoherence of the spins was studied using the kinetic 
equation for the reduced density matrix. Additionally, 
a quantum gate consisting of two quantum dots embed- 
ded in a two-dimensional electron gas of GaAs/AlGaAs 
heterostructure was investigated. The RKKY effect was 
provided through the s — d exchange interaction which 
was assumed to take place just at the positions of the 
localized spins, that is the interaction radius was zero. 
Tamura et al~2. studied the RKKY interaction between 
the localized spins of two quantum dots placed at the op- 
posite edges of a one-dimensional (ID) conducting chan- 
nel. The RKKY interaction between the spins across the 
channel was taken into account by virtue of an exchange 
interaction where the Fourier transform of the coupling 
constant was assumed to be momentum-dependent. This 
in principle means that the interaction could be non- 
local. However, consequences of this non-locality were 
out of focus of that work. In Ref. [H the RKKY inter- 
action in a coupled quantum dot system embedded in 
a ring with a spin-orbit interaction was explored in the 
presence of the Aharonov-Bohm and Aharonov-Casher 
effects. The s — d exchange interaction responsible for 
the RKKY effect in this system was also local. 

Finally, we would like to mention that effects of the 
finite size of the ionic spin distribution on the RKKY in- 
teraction have been studied in bulk systems in connection 
with ferromagnetic Heusler alloy o 14 i 15 ' 16 i 17 ' 18 . However, 
the models used in those attempts did not give the delta 
function in the limit ro — > 0. 

The purpose of the present work is to verify by means 
of a simple model whether in a mesoscopic ring a finite 
value of ro can play any role in the formation of the elec- 
tron spin density oscillations which underpin the RKKY 
indirect exchange interaction. In the limit ro — > the 
model which we use gives the delta function. It is shown 
that if the interaction radius is finite, its impact on the 



electron spin density oscillations is to reduce the oscilla- 
tion amplitude. An interesting feature is that when the 
number of the electrons in the gas is large enough, this 
reduction, produced by the presence of a small domain of 
size ro , is significant and takes place in the whole system 
with much larger size. 

The paper is organized as follows. Sec. [Til describes a 
mathematical formulation of the problem. In Sec. Illll thc 
solution is obtained using the Matsubara diagrammatic 
approach. The results are discussed in Sec. IIVI 



II. THEORETICAL MODEL 

We consider a non-interacting Fermi-gas on a ID ring 
of radius R. The electron positions are specified by the 
polar angle </>. An ion (or another object) with a non-zero 
total spin is placed on the ring at <j) = 0. The system is 
schematically shown in Fig. [1] The interaction between 
the localized spin and conduction electrons is assumed 
to be an exchange interaction which in general can be 
non-local. 

To formulate the problem mathematically we write 
down the Hamiltonian of the system in the form: 



H = H n 



inl • 



(i) 



In Eq. (UJ) Ho is the Hamiltonian of the non-interacting 
TV-electron system on a ID ring: 



Ho 
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with 



H ° 2mR 2 ' 



(2) 



(3) 



where m is the electron mass and fa is the- i-th elec- 
tron coordinate. The spin-degenerate eigen-values of the 
single-particle Hamiltonian Hq are 



(4) 



2mR 2 ' 



where n = 0, ±1, . . . and a is the spin index. The nor- 
malized eigen-states of Hq, \na), in the coordinate rep- 
resentation are 



(<j)a'\na) = S a 



1 



exp( 



(5) 



The term H lnt in Eq. (Q} describes the exchange interac- 
tion between the localized spin and the electron gas and 
it is conventionally written as: 



N 

E 

3=1 



J(fa)S l a), 



(6) 
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FIG. 3: (Color online) The relative change of the spin density 
at 4> = 7T (the furthermost point from the localized spin). The 
reference is the spin density at ro = 0. The temperature is 
T = IK. 



the ring is defined as 



N 



(10) 



The electron spin density at point <j) is obtained through 
the statistical average using the Gibbs grand canonical 
ensemble: 



K(<t>, r o) = 



Tr [cxp ( 



H-nN \ 



(ii) 



where we explicitly show that the electron spin density is 
also a function of the interaction radius tq. In Eq. pip 
Z is the partition function: 



Z = Tr 



exp 



H- fiN 
knT 



(12) 



where J(<f>) is the coupling function of the exchange in- 



teraction, S is the localized spin placed at 



0, a, is 



the vector of the electron spin Pauli matrices and the 
summation over the index i is assumed. 

As it was mentioned above, traditionally it is assumed 
that the exchange interaction is local, that is the embed- 
ded spin interacts with the electrons only at <j> = 0. In 
this case the coupling function </(</>) is modeled by the 
following dependence on the polar angle: 



J(#) = J<5(0), 



(7) 



where J is a coupling constant. 

In this work we suggest a simple model in which the lo- 
calized spin interacts with the electrons in a small vicinity 
of the point = and when the vicinity is tightened into 
a point at = 0, the model turns into the conventional 
one ©: 



J(0) 




'(to) ]' 



'2tt- 



< 7T, 

7r < 6 < 2tt. 



(8) 



The size 0o = ro/R of this vicinity is estimated from 
the radius ro of the ion (or from a characteristic size of 
another object) which is the source of the spin centered 
around = 0. Using the well known representation of 
Dirac's delta function, 



S(x) = lim 
t-»o 



-p[-(7?) 2 ] 



Eq. 



(9) 



takes 



one is convinced that the non-local model, 
the local form, Eq. jJJ, when r — > 0. 

To study whether finite values of ro have any effect we 
will calculate the electron spin density at finite temper- 
atures. The electron spin density operator at point on 



In Eqs. (jTTJ) and (fT!?|) /c B is the Boltzmann constant, \i 
is the chemical potential, T is the temperature, N is the 
particle number operator and the trace is taken using 
a complete set of states of the Fock space. It is also 
important to note that since the number TV of electrons 
is fixed, the chemical potential fi is not an independent 
variable but a function of the temperature T. 



III. IMAGINARY TIME GREEN'S FUNCTION 
SOLUTION 



The problem formulated in the previous section is ob- 
viously not solvable exactly. Thus some approximation 
methods should be applied. For small values of the cou- 
pling constant J in Eq. ([5]) a perturbation theory, where 
.Hint is considered as a perturbation, can be applied. As 
the calculations are performed at finite temperatures, in- 
stead of the pure quantum mechanical perturbation the- 
ory one has to use the so-called thermodynamic perturba- 
tion theory^ for the quantum statistical averages. This 
perturbation theory is in general valid when the pertur- 
bation energy per particle is less than k^T, i.e., J < kBT. 
However, very often it happens that when T — > 0, the co- 
efficients of the perturbation expansion change as func- 
tions of T in such a way that the thermodynamic pertur- 
bation theory can remain valid for all temperatures. 

Although the thermodynamic perturbation theory 
gives a general approach to calculate statistical averages, 
in its original form it is quite cumbersome. It is more 
convenient to use this theory reformulated in terms of a 
diagrammatic approach, namely the Matsubara (or imag- 
inary time) diagrammatic method^. 

In order to employ this technique for our purposes we 
first rewrite the total Hamiltonian of the problem, Eq. 
([T|), in the second quantized form using the {na} single- 
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particle basis: 
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x exp[j0(n - n')] \al U7 a n > a 



(13) 



The second quantized form of the electron spin density 
operator at point is 

K{<t>) = ir- (^V)exp[#(n-n')]aL<W (14) 

Z7T * — ' 

na n' a' 

and the expression for the electron spin density at point 
may be rewritten as 



n s {(p,r 



)=-^2{a'\a i \a)g aa ,(4>,T;4>,r). (15) 



In the last expression Q aa i ((j), r; 0' , r') is the one-particle 
imaginary time Green's function defined as 



<w(0,t;0',t') 



;Tr 



exp 



k*T 



(16) 



where T is the time-ordering operator and CT (0, r) are 
the imaginary time field operators, 



(17) 



0v(0,t) = 

= exp[r(# - nft)]$„(4>) exp[-r(# - uN)}, 
with ip/rift) related to the annihilation operators a na as 
ipa{(j>) = -^=y^exp(-z0n)a TI(T . (18) 



2tt 

IL 

We now apply the diagrammatic expansion of the Green's 
function Gcrcr'(<p, t;0',t / ). The effect of the RKKY spin 
density oscillations appears already in the first order and 
thus we only need to consider the first order diagram. 
Such an approach to the RKKY spin density oscillations 
was considered in Ref. l2l| for a three-dimensional elec- 
tron gas. However, in that case the electron spectrum 
was continuous and to perform momentum integrals the 
linearization of the spectrum on the Fermi-surface was 
employed to get the long range behavior of the RKKY 
oscillations. In our problem the electron spectrum is dis- 
crete and instead of integrals we will have sums over dis- 
crete indices. Moreover, since our system is finite we are 
interested in the RKKY oscillations in the whole range 
and not only at long distances from the localized spin. 
The first order contribution to the Green's function 

Gaa'(<f>,T;<f>',T) is 



^(0,t;0',t) 



JS* 
(2ir) 2 H { 



x ) exp(—iri(f) + in' ft ) jn n' $n n' ? 



(19) 



n n 
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FIG. 4: (Color online) The electron spin density oscillations 
in the vicinity of the point cj> = 7T. The number of electrons is 
Nei = 600. The temperature is T = IK. 



where 



and 



Jn n' — exp 

x Re I erf 



-(n-n') 2 4, 2 



7T I . 

— + -(n- 7i)0o 

00 2 



(20) 



n 12 ± n 2 , 

i>nn' — \ ™ , " k ,n n 

[-n„(l - rin)-^, n' 2 = n z 

with n n being the fermion occupation numbers, 

1 



exp(^) + 1 



Since 4>q <C 7T, we have 

Jn n' ~ exp 



--An-n') 2 ^ 



(21) 



(22) 



(23) 



and this approximation will be used in the calculations 
below. 

Substituting Eq. (|19[) into Eq. (| 1 5|) and taking into 
account Eq. (|2ip . we obtain the first order contribution 
to the electron spin density 



ni(4>,r ) 



2JS* 



n B (<f>,r ), 



(24) 



where n s (<f>, rn) is given by the following expression: 
n s (0,r o ) = ^2 exp[-i(n- n')<j>]j nn > _ " - 

n n 
(«'V» 2 ) 

- ^2 exp[-i(n-ri)<f>]j n . 

n n 
(n' 2 =n 2 ) 



n n (l - n n ) 



(25) 
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For a given value of the interaction radius ro the function 
n s{ ( t ) i r o) provides an oscillatory behavior of the electron 
spin density as a function of the polar angle <f>. 



IV. RESULTS AND DISCUSSION 

In this section we numerically analyze the electron spin 
density oscillations. Using Eq. ([25]) we calculate the 
function n s (<j), ro). The mesoscopic ring is assumed to be 
fabricated on AlGaAs-GaAs heterostructures. The val- 
ues of the parameters are taken close to the ones used in 
previous works^ and in experiments^. The ring radius 
is R = 40 nm, the effective mass m = 0.067m , where 
mo is the free electron mass. 

Let us first consider the behavior of n s {4>, ro) as a func- 
tion of the polar angle for the conventional model with 
ro = 0. It is shown in Fig. [5] for different values of the 
temperature T and for the number of electrons N c \ = 24. 
In systems with a continuous spectrum the RKKY spin 
density oscillations behave like cos(2ppr/fi,) (similar to 
the Friedel oscillations of the electronic density) where 
Pp is the Fermi momentum and r is the distance from 
the localized spin. In our case the spectrum is discrete. 
The analog of the Fermi momentum is the energy level 
number np such that at T = the states with \n\ > np 
are not populated. Since the state with a given value of 
n can be occupied by two electrons, for iV e i = 24 one gets 
np = 6 and thus 2np = 12. In agreement with this esti- 
mation Fig. [2] shows 12 oscillations. As it was discussed 
in the literature (see, for example, Ref . l2ll) . the effect of 
the temperature on the RKKY spin density oscillations is 
to produce a faster reduction of the oscillation amplitude 
with the distance from the localized spin. For example in 
three-dimensional electron gases at T = the oscillation 
amplitude at large distances decreases as 1/r 3 and at fi- 
nite temperatures it is damped at the thermal distance 
Hpp I \2Timkp,T) . An analogous behavior takes place in 
our case as well. As it can be seen from Fig. [2] the ampli- 
tude of the oscillations as a function of the polar angle is 
damped faster for higher values of the temperature. 

Next we turn our attention to the effect of the finite 
interaction radius ro of the exchange interaction on the 
electron spin density oscillations. From Eq. (|23|) it seems 
that for realistic, that is small, values of <f>o the function 
jnn' ~ 1 and a finite interaction radius value does not 
produce any change in the electron spin density oscil- 
lations in comparison with the case r = 0. However, 
this reasoning is not entirely correct because it does not 
take into consideration the number of electrons in the 
gas surrounding the localized spin in the ring. Indeed, 
when the number of electrons in the gas grows, energy 
levels with higher values of \n\ become important. The 
contributions with higher values of \n — n'\ are involved. 
From Eqs. (|23[) and ([25]) one observes that in parallel 
with this involvement the contributions from terms with 
high values of \n — n'\ get more suppressed. It does not 
matter how small the interaction radius is. The main 



point is that it is finite. For any finite value of ro there 
exists a number of electrons N c in the gas such that for 
iV c i > A c the interaction radius, no matter how small it 
is, will produce more and more pronounced impact on 
the spin density. This is demonstrated in Fig. [3] where 
the relative change of the spin density at (j> = tt is de- 
picted as a function of the number of the electrons in 
the gas for different values of the interaction radius ro- 
As expected the change of the electron spin density in 
comparison with the conventional model with r = is 
negligible for small numbers of the electron in the ring. 
For larger numbers of the electrons the electron spin den- 
sity for finite values of ro starts to deviate from the model 
with ro = 0. It is interesting that even the interaction 
radius ro = 0.05 nm can produce an observable change 
of approximately 11% for N c \ — 1600. The ring distance 
between the two points = and <p = n is about 126 nm 
while the size of the area where the spin is localized is 0.1 
nm, that is three orders of magnitude less. An important 
result is that the small vicinity around the localized spin 
is able to significantly change the electron spin density in 
every point of the system whose size is several orders of 
magnitude larger than the size of the domain over which 
the localized spin is spread. 

The RKKY spin density oscillations for N c \ = 600 are 
displayed in Fig. |4] for different values of ro. In this case 
the number of the oscillations is approximately equal to 
300 and thus only a small vicinity around <f> = tt is plotted 
to clearly show the oscillations. 

Finally, we note that semi-quantitatively the oscillat- 
ing behavior can be explained by the dominance of the 
term cos(2n F 0)j„ F: _„ F n„ F (1 - n„ F ) in Eq. (gSJ). The 
2np RKKY oscillations with ro = are weighted with 
jn F ,-n F < 1. For small N c \ the weight j UF ,-n F ~ 1 and 
plays no role but for large N e \ this weight reduces the 
amplitude of the 2np RKKY oscillations by the factor 

CXp[-(?lF<M 2 ]. 



V. SUMMARY 

The Ruderman-Kittel-Kasuya-Yosida (RKKY) spin 
density oscillations, induced by an exchange interaction 
between a localized spin and the electron gas in which 
the spin is embedded, have been investigated at finite 
temperatures taking into account finite values of the ex- 
change interaction radius. It has been found that the 
electron spin density in a non-interacting gas on a meso- 
scopic ring is changed in comparison with the traditional 
model which assumes that the interaction radius is zero. 
The amplitude of the RKKY oscillations is suppressed 
when the number of the electrons in the gas increases. 
This suppression gets stronger for larger values of the 
interaction radius. 

A remarkable point is that as soon as the interaction 
radius is finite, it is already not important how small it is 
because its impact always becomes observable when the 
number of the electrons in the ring is large enough. 
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